LPN 11-91 



Z"-boson production in association with a top anti-top pair at 
NLO accuracy with part on shower effects 

M. V. Garzelli/'^ A. Kardos/'^ C. G. Papadopoulos,^ and Z. Trocsanyi^'^ 

^Institute of Physics, University of Debrecen, 
H-4010 Debrecen P.O.Box 105, Hungary 
"^Laboratory for Astroparticle Physics, 
University of Nova Gorica, SI- 5000 Nova Gorica, Slovenia 
^Institute of Nuclear Research of the Hungarian Academy of Sciences, Hungary 
^Institute of Nuclear Physics, NCSR Demokritos GR-15310 Athens, Greece 

(Dated: December 22, 2011) 

Abstract 

We present predictions for the production cross section of a Standard Model Z*^-boson in associa- 
tion with a tt pair at the next-to-leading order accuracy in QCD, matched with shower Monte Carlo 
programs to evolve the system down to the hadronization energy scale. We adopt a framework 
based on three well established numerical codes, namely the POWHEG-BOX, used for computing the 
cross section, HELAC-NLO, which generates all necessary input matrix elements, and finally a par- 
ton shower program, such as PYTHIA or HERWIG, which allows for including t-quark and Z^-boson 
decays at the leading order accuracy and generates shower emissions, hadronization and hadron 
decays. 



PACS numbers: 12.38.-t, 13.87.-a, 14.65.Ha, 14.80.Hp 



With increasing collider energies, the t-quark plays an increasingly important role in 
particle physics. Its production cross section grows faster with energy than that of any 
other discovered Standard Model (SM) particle. Already after the first year of successful 
run of the LHC, the tt production cross section is measured with unprecedented accuracy 
at ^/s = 7TeV, so that the corresponding SM theoretical prediction will be challenged 
soon [H [2] . However, many other t-quark properties have not yet been directly accessed. In 
particular, its couplings to neutral gauge (especially the and scalar bosons are still prone 
to large uncertainties. In Refs. [31 E] the possibility of measuring the ttZ and tt7 couplings 
was studied based upon leading-order (LO) parton level predictions. Although such precision 
is sufficient for feasibility studies, finding the optimal values of the experimental cuts requires 
indeed predictions at higher accuracy. 

An essential step towards higher accuracy is the inclusion of next-to-leading order (NLO) 
radiative corrections. Recent theoretical advances made possible our computation of the 
pp — 7- ttZ cross section at the parton level, including QCD corrections at NLO [5j. In 
order though to get the optimum benefit and to produce predictions that can be directly 
compared to experimental data at the hadron level, a matching with parton shower (PS) and 
hadronization implemented in shower Monte Carlo (SMC) programs is ultimately inevitable. 
Thus, in this letter we present first predictions for pp ttZ production at LHC at the 
matched NLO + PS accuracy. 

In constructing a general interface of PS to matrix element (ME) computations with 
NLO accuracy in QCD, we have chosen to combine the POWHEG |7| method and FKS 
subtraction scheme [8], as implemented in the POWHEG-BOX [9j computer framework, with the 
HELAC-NLO [10] approach, respectively. In particular, POWHEG-BOX requires the relevant 
MEs as external input. We obtain the latter in a semi-automatic way by codes in the 
HELAC-NLO package [11] . With this input POWHEG-BOX is used to generate events at the Born 
plus first radiation emission level, stored in Les Houches Event Files (LHEF) [12], that can 
be interfaced to standard SMC programs. Previous applications of the whole framework, 
proving its robustness, were presented in Refs. [T3l E]. This same setup also allows for 
exact NLO pure hard-scattering predictions. Further details on the implementation of the 
computation of the pp — i- ttZ hard-scattering cross-section in it, at NLO accuracy in QCD, 
together with checks, were recorded in Ref. [S]. 
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Figure 1. Transverse momentum (left) and rapidity (right) distributions of the Z^-boson at NLO 
and after first radiation (PowHel). The lower panels show the ratio of the two predictions with 
combined uncertainties. 



All these computations are steps of an ongoing project for generating event samples for 
pp — 7- ttX processes, where X stays for a hard partonic object. The events we generate 
are stored in LHEF, made available on the web, and are ready to be interfaced to standard 
SMC programs to produce predictions for distributions at the hadron level. Such predic- 
tions can be useful for optimizing the selection cuts applied to disentangle the signal from 
the background, in order to improve the experimental accuracy of the t-quark coupling 
measurements. 

Interfacing NLO calculations to SMC programs allows to estimate the effects of decays, 
shower emissions and hadronization, therefore we have analyzed the process at hand at three 
different stages of evolution: 

PowHel: we analyzed the events including no more parton emissions than the first 
and hardest one, collected in LHEF produced as output of POWHEG-BOX+ HELAC-NLO 
(PowHel). 

Decay: we just included on-shell decays of t-quarks and the Z°-boson, as implemented 
in PYTHIA [15J, and further decays of their decay products, like charged leptons (the r 
is considered as unstable) and gauge bosons (W), turning off any initial and final PS 
and hadronization effect. 
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Full SMC: decays, shower evolution, hadronization and hadron decays have been 
included in our simulations, using both PYTHIA and HERWIG [TB] . 

In our computation, we adopted the following parameters: ^/s = 7TeV, CTEQ6.6M 
PDF set from LHAPDF, with a 2-loop running a^, 5 light flavours and A^^ = 226 MeV, 
nit = 172.9 GeV, mz = 91.1876 GeV, Gp = 1.16639 ■ 10"^ GeV~^ The renormalization 
and factorization scales were chosen equal to the default scale = rrit + mz/l. We used 
the last version of the SMC fortran codes: PYTHIA 6.425 and HERWIG 6.520. Following our 
implementation of \XH hadroproduction in Ref. [11], in both SMC setup muons (default in 
PYTHIA) and neutral pions were assumed as stable particles. All other particles and hadrons 
were allowed to be stable or to decay according to the default implementation of each SMC. 
Masses and total decay widths of the elementary particles were tuned to the same values in 
PYTHIA and HERWIG, but each of the two codes was allowed to compute autonomously partial 
branching fractions in different decay channels for all unstable particles and hadrons. Mul- 
tiparticle interaction effects were neglected (default in HERWIG). Additionally, the intrinsic 
p_L-spreading of valence partons in incoming hadrons in HERWIG was assumed to be 2.5 GeV. 

First, to check event generation, we compared several distributions from events including 
no more than first radiation emission (PowHel level) with the NLO predictions of Ref. [3]. 
We found agreement for all considered distributions. As examples, we show in Fig.[T] the 
transverse momentum and rapidity distributions of the Z°-boson. 

Next, we studied the effect of the full SMC by comparing distributions at the decay and 
SMC level. Since particle yields are very different at the end of these two stages, we made 
such a comparison without any selection cut, in order to avoid the introduction of any bias. 
As an illustrative example, we present the distributions of the transverse momentum and 
rapidity of the hardest jet, p-^ and Uj, in Fig.[2| Jets are reconstructed through the anti — k± 
algorithm with R = 0.4, as implemented in FastJet [1^. The softening of the transverse 
momentum spectrum is apparent as going from the decay level to the full SMC one, while 
the effect of the shower on the rapidity of the hardest jet is almost negligible and rather 
homogeneous. The cross-section at both level amounts to cr = 138.7 ± 0.01 fb. Using our 
setup for the full SMC's, we found agreement between PYTHIA and HERWIG predictions within 
very few percent, despite the conceptual differences between the two SMC generators as for 
the shower ordering variables and hadronization models, confirming the level of agreement 
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Figure 2. Transverse momentum (left) and rapidity (right) distributions of the hardest jet after 
decay and after full SMC. The lower panels show the ratio of all predictions to PowHel+SMC using 
PYTHIA. 



already reported in Ref. [H] in the study of a different process. 

Next, we made predictions for ttZ hadroproduction at the LHC including experimental 
selection cuts. For this analysis, in the absence of a dedicated tune for NLO matched 
computations, PYTHIA was tuned to the Perugia 2011 set of values, one of the most recent 
LO tunes jTS], updated on the basis of recent LHC data, providing a p^-ordered PS. Its 
application turned out to increase our particle yields by about 10%. As a consequence, 
the agreement between the tuned PYTHIA and untuned HERWIG predictions decreases (as for 
HERWIG, the default configuration was used, providing instead an angular-ordered PS), and 
we present only the PYTHIA ones. 

In case of ttZ hadroproduction overwhelming backgrounds come from tt+jets final states. 
In Ref. |1] the differential cross section as a function of missing transverse momentum for 
the production of p±bb+4: jets was found a useful tool for differentiating the signal and the 
possible backgrounds. The proposed set of selection cuts is rather exclusive and the rates 
decrease so much that the measurement for the present LHC run at ^/s = 7TeV looks 
quite demanding from the statistical point of view, therefore, we restrict ourself to present 
predictions for the future runs at y/s = 14TeV (craiicut,i4/'7aiicut,7 ~ 7 and 8 at the decay 
and at the full SMC level, respectively). 

In Fig. [3] we show the distributions of transverse momentum and rapidity of the hardest 
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Figure 3. Transverse momentum (left) and rapidity (right) distributions of the hardest jet after 
decay and after full SMC (PYTHIA), under selection cuts (1-8) implemented at both levels. The 
lower panels show the ratio of the predictions at different levels. 



jet using the following reduced set of cuts: 1) we reconstruct at least six jets with rapidity 
\y\ < 2.5, 2) of these we require at least one 6-jet and one 6-jet, 3) for 6-jets p5_ > 20 GeV, 
4) for other jets p^""" ^ > 30 GeV, 5) at least 3 jets {b or non-6) with p{ > 50 GeV, 6) 
AR(j,j) > 0.4, where j denotes any (6 or non-6) jet and AR is defined as a/ A</)^ + Ay^ , 
7-8) A0(|i^,p_Lj) > 100°, with meaning either (p_l(6i) +^±(62)) (cut 7), or (p_L(ji) + 
P±(j2) +P±(j3) +P±(j4)) (cut 8), where 61, 62 and ji, j2, Js, J4 are the jets that allow for 
the best t — > bW+ — t- hjj and t— )■ hW~ — )■ hjj invariant mass simultaneous reconstruction, 
since they minimize the 

^2., ■ . 7 . ^. ^ _ (^jii2 - ^w? (m,3,, - mwf 
X [O1J132, O2J3J4) — 2 ^ 2 

computed by considering all possible jfcji, 6jjfcji and hijkji combinations. The W — )■ jj 
and t— 7- 6jj invariant mass resolutions were set to cxvy = 7.8 GeV and Ot = 13.4 GeV , 
respectively [19J. The PowHel+PYTHIA cross sections after these cuts amount to a^ec = 
65.56 ±0. 15 fb and ctsmc = 53.74 ± 0.13 fb. 

In Fig.|4]we plot the invariant mass distribution of the t-quark, as reconstructed from its 
decay products, by minimizing the above. At the decay level, the reconstruction leads to 
a clear peak centered around the rrit value. On the other hand, after full SMC, due both to 
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Figure 4. Invariant mass distribution of the t-quark reconstructed from the decay products at both 
decay and full SMC levels, for the ttZ signal and, at the decay level, for one background (tt+jet) 
after selection cuts (1-8) (wider distributions in abscissa values) and after selection cuts (1~10) 
(narrower distributions). 



further emissions which modify jet content and to hadron decays, there are more candidate 
jets and the reconstruction is less successful. Although a peak is still visible (more evident 
in non-log scale), it is smeared towards lower mass values. The effect of the shower and 
hadronization turns out to be especially large in the peak region. 

In Fig. [4] we also show the m^jj distribution after decay for an important background 
process: tt-pair production associated with a jet (obtained at the scale /io = rrit). Clearly, 
the background overwhelms the signal, therefore, in order to select the peak region, we 
include two more cuts: 9) j[)±{due to all z/'s) > (of J^ts, b or non-6), and 10) 

Xmin < 3, where Xmin is the minimum of the above. Thus, we closely reproduce the cuts 
in Ref. [1], aimed at favoring the — )■ uu decay channel. The effect of the whole set of 
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Figure 5. Distribution of the missing transverse momentum after decay, under physical cuts (1~10) 
applied to the signal (tiZ) and to one background (tt+jet). 



cuts on top reconstruction in tiZ and tt+jet events is also shown in Fig.|4| Although this 
set of cuts is effective in selecting the signal, the background is globally still larger: for the 
signal (Jdec = 4.83 ± 0.04 fb, while for the background adec = 9.86 ± 1.05 fb. However, as 
can be understood from Fig.|5} where the distributions of the missing transverse momentum 
after decay are shown for both ttZ and tt+jet, these cuts allow for disentangling the signal, 
at least at the decay level. At the shower level, the f)± distributions of the ttZ signal still 
shows a harder spectrum than the one of the tt+jet background, but to a lesser extent. In 
this case, the effect of different top reconstruction strategies, still under investigation, can 
be crucial to help better disentangle the signal from the background in the ^_l66+4 jets 
considered channel. 

We studied the hadroproduction of a Z° boson in association with a tt-pair, process of 
interest for measuring the ttZ-couphng directly at the LHC. We studied the effect of heavy 
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particle decays as well as the one of the full SMC. We produced predictions for the LHC. 
As the production cross section is rather small, measuring the ttZ-coupling becomes more 
feasible after the planned 14TeV energy upgrade. Once all background processes will be 
predicted with the same accuracy, our predictions will make possible a realistic optimization 
of the experimental cuts. 
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